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Abstract 

The  fundamental  solution  for  VLF  waves  propagating 
past  a  shoreline  In  the  earth-Ionosphere  duct  Is  obtained 
In  terms  of  an  Infinite  series  using  the  Wlener-Hopf  methods. 
The  effects  of  the  land,  sea  and  Ionosphere  are  given  In 
terms  of  complex  Impedance  boundary  conditions.   As  a  result, 
each  mode  in  the  series  is  progressively  damped.   At  low 
frequencies  and  distances  far  from  the  shoreline  where  a 
single  mode  dominates,  reflection  and  transmission  coeffi- 
cients are  obtained.   When  the  ionospheric  Impedance  is 
small  and  when  the  source  and  receiver  are  sufficiently 
elevated,  the  effect  of  the  shoreline  is  not  observable. 
However,  when  either  the  source  or  receiver  is  near  the 
earth  the  portion  of  the  wave  reflected  is  proportional  to 
the  difference  of  sea  and  land  Impedances.   At  slightly 
higher  frequencies,  there  exists  an  annular  region  about 
the  shoreline  where  two  modes  propagate  over  the  sea  while, 
due  to  heavy  damping,  only  one  mode  propagates  over  the 
land.   The  resulting  increase  of  the  field  strength  as  the 
receiver  passes  from  land,  seaward  is  the  "recovery  effect." 
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1 .   Introduction 

The  purpose  of  this  paper  is  to  derive  and  analyze 
the  fundamental  solution  (Green's  function)  for  very-low- 
frequency  (VLP)  two-dimensional  electromagnetic  waves 
propagating  past  a  shoreline  in  the  earth-ionosphere  wave 
guide.   Although  portions  of  this  paper  have  appeared  as 
a  report  (Dobrott  1960),  continued  interest  in  the  subject 
(Wait  1965)  suggests  that  these  results  be  made  more  readily 
available. 

The  original  version,  which  gives  the  fundamental 
solution  in  terms  of  an  infinite  set  of  spatiaJly  damped 
modes,  has  been  condensed.   New  results  are  given  on  the 
reflection  and  transmission  coefficients,  when  on]y  one  of 
these  modes  propagates,  and  on  the  "recovery  effect" 
(Mlllington  19^9)?  where  more  than  one  mode  propagates 
resulting  in  an  increase  in  the  magnitude  of  the  vertical 
component  of  the  electric  field  as  the  wave  passes  from 
land  to  sea.   Wiener-Hopf  techniques  are  used  to  solve  the 
appropriate  partial  differential  equations  with  mixed 
boundary  conditions.   In  particular,  this  paper  follows 
methods  due  to  D.S.  Jones  as  presented  by  Noble  (]9'''8). 

The  analysis  of  electromagnetic  wave  propagation  over 
land-sea  boundaries  has  an  old  and  extensive  literature. 
These  papers,  including  some  which  utilize  Wiener-Hopf 
methods,  are  not  generally  applicable  to  VLP  propagation 
because  of  the  absence  of  the  ionospheric  boundary  in  their 
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analysis.   A  partial  list  of  this  Indirectly  related  liter- 
ature Is  given  at  the  end  of  this  paper. 

Several  papers  using  the  Wlener-Hopf  technique  but  not 
dealing  with  electromagnetic  wave  propagation,  exhibit 
mathematical  similarities  to  VLP  mode  transmission.   Heins 
and  Peshbach  (19^7)  studied  the  coupling  of  sound  waves 
between  two  dissimilar  acoustical  ducts.   Helns  (19^8) 
extended  these  results  when  studying  the  behavior  of  water 
waves  In  a  channel  of  finite  depth  Incident  upon  a  dock. 
Similarly,  Weitz  and  Keller  (1950,1953)  have  examined  the 
reflection  of  water  waves  from  floating  ice  in  water  of 
finite  depth.   The  present  paper  differs  from  the  above  in 
that  all  three  boundaries,  i.e.  ionosphere,  land  and  sea, 
are  represented  by  complex  impedances.   In  the  above  papers, 
the  equivalent  of  the  boundary  impedance  corresponding  to 
the  ionosphere  and  sometimes  that  of  one  other  boundary  has 
been  set  equal  to  zero.   Such  a  restriction  on  the  ionospheric 
boundary  is  physically  unreasonable  since  the  nature  of  the 
waves  in  the  earth-ionosphere  duct  is  strongly  determined 
by  all  the  boundaries. 

Figure  1  is  a  schematic  of  the  problem  to  be  considered. 
The  z-axis  corresponds  to  the  shore  line  in  a  right-hand 
rectangular  coordinate  system.  The  earth  surface  lies  in 
the  x-z  plane  and  the  ionosphere  is  taken  to  have  a  sharp 
boundary  at  y  =  0.  The  earth's  magnetic  field  is  assumed 
parallel  to  the  shore  line.  The  symbol  T  indicates  the 
source  and  R  is  the  receiver  or  field  point.   The  sources  and 


the  fields  are  assumed  to  vary  In  time  like  exp(lcot). 
Appropriate  to  the  reduction  to  a  two-dimensional  problem, 
the  source  position  and  strength  is  assumed  independent  of 
the  z-coordinate. 

The  surface  impedance  representing  the  effect  of  a 
sharply  bounded  ionosphere  depends  upon  the  direction  of 
propagation  relative  to  the  orientation  of  the  earth's 
magnetic  field.   For  propagation  across  the  earth's  field. 
Maxwell's  equations  rigorously  decouple  into  transverse 
magnetic  (TM)  and  transverse  electric  (TE)  waves.   For  low 
frequencies,  the  ionospheric  impedance  is  given  approximately 
by 


2 
cu  +v 

2 

H 

exp 


TT 


„    tan   (— ) 
2         ^Od^ 


(i.l) 


where 


_1 X 

Z   H 
o   z 


_1 z 

Z^  H 
o   X 


for  TM  waves 


y  =  o 


y  =  o 


for  TE  waves 


(1.2) 


where  Z^  =  377  ohms  is  the  characteristic  impedance  of  free 
space,  OJ  the  plasma  frequency  and  v  the  electron  collision 
frequency.   A  signal  frequency  cu  is  considered  low  when 
CD  «  GO  «  CO  ,  where  co   is  the  electron  cyclotron  frequency 
in  the  earth's  field.   At  higher  frequencies  the  ionospheric 
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impedance  for  TE  and  TM  waves  are  different  (see  Dobrott 
1961  )j  however  the  fundamental,  solution  may  be  trivially- 
modified  by  substitution  of  the  appropriate  expression  for 
Z^  in  the  separate  cases. 

The  relative  land  and  sea  Impedances  are  given  by 
(see  Kerr  1951) 


'L,S 


+ 


where 


^L,S 


J 


■L,S 


_1 X 

Z   H 
o   z 


exp 


1   ^L,S 

'  ^^0^1 


L,S 


(1.3) 


for  TM  waves 


y  =  o 


(1.4) 


^0  ^x 


for  TE  waves 


y  =  o 

where  e-,     is  the  relative  dielectric  constant  and  a^ 
is  the  conductivity  for  land  and  sea  respectively.   All 
impedance  magnitudes  are  assumed  small,  in  the  order 
IZ3I  «  |Zj|  «  |Z^|  «  1. 

The  z  components  of  the  electric  and  magnetic  fields 
produced  by  a  current  distribution  j(x,y)  exp(icot)  are 
given  in  terms  of  the  fundamental  solution  g(x  ,y    |x,y)  as 

H^(x,y,t)  =  exp(i(Dt)  <^     J_(xQ,y^  )xV^.  i^g(xQ,y^|  x,y )  dx^dy^ 


source 


where  1   is  the  unit  vector  in  the  z-direction  and 
z 


(1-5) 
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E^(x,y,t)  =-ia)^^exp(icut 


)  §  i(=<o- 


V  )•!  efx  ,y  |x,y)dx  dy 


source 


(1.6) 


where  the  fundamental  solution  g(x^,y^  x,y)  satisfies  the 


equation 


(1.7) 


for  all  X  4=  X   and  y  +  y^,  the  condition 


lim  ^ 


f 

^0+^ 

r- 

dy 

"Sx 

^0-^ 

„ 

V 

X  +e 
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Sg 


^o-^ 
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dx 


X  -e 
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^g 


^ 


y  +e 
''o 


Sg 


and  the  boundary  conditions 


4^  +  ikZ^g  =  0   when   y  =  b. 


<  X  < 


^o-^ 


^ 


(1.8) 


(1.9) 


1^  -  IkZ  g  =  0   when   y  =  0,  0  <  x  < 


(1.10) 


ll  -  IkZgg  =  0   when   y  =  0,  -c»  <  x  <  0. 


(1.11) 


where  k  =  oo/c  and  c  is  the  velocity  of  light  in  a  vacuum. 

For  the  purposes  of  uniqueness,  it  is  necessary  to 
specify  conditions  at  the  shoreline  and  at  infinity.   At 
the  shoreline,  electromagnetic  power  is  neither  absorbed 
nor  radiated.   This  is  equivalent  to  the  statement. 
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11m 
p->q 


<~  ^                                                       ag(x  ,y  Ip  cos  0,p  sin  0) 
g  (x^,yQ|p  cos  0,p  sin  0)  '± — "^ ^- pd0  =  0 


(1.12) 

* 

where  g  denotes  the  complex  conjugate  of  g,  x  =  p  cos  Q, 

y  «  p  sin  Q    and  0  <  0  <  it.   Similarly,  since  the  duct 
boundaries  are  lossy,  no  power  will  reach  infinity.   This 
statement  is  the  same  as  saying 


lim 
x|— ».<» 


b 
''  *.  ,    ,  Sg(xQ,yQ|x,y) 


g  (x^,y^|x,y)  "^-^ dy  =  0.  (1-13) 


o 

the  z  components  suffice  to  characterize  the  TM  and 
TE  modes  for  this  problem.   The  other  components  are 
determined,  in  principle,  by  Maxwell's  equations.   Physically, 
ia  g(x  ,y  |x,y)exp(lcot)  is  the  two-dimensional  vector 
potential  at  (x,y)  produced  by  a  dlpole  source  at  (x  ,y  ). 

In  Section  2,  the  integral  representation  for  the 
fundamental  solution  is  given  in  terms  of  an  integral 
involving  the  Fourier  transform  in  x  of  g  evaluated  at  y  =  0. 
This  transformed  function  is  unknown  and  is  determined  in  an 
auxiliary  problem  by  Wlener-Hopf  techniques  (Section  3).   An 
explicit  series  representation  for  the  fundamental  solution 
is  given  in  Section  4.   The  particular  asymptotic  form  for 
the  solution  well  away  from  the  shoreline  is  used  to  obtain 
reflection  and  transmission  coefficients.   The  recovery 
effect  is  obtained  by  examining  the  solution  in  an  annulus 
about  the  earth-sea  juncture  for  frequencies  which  permit  more 
than  one  mode  to  propagate  in  the  earth-ionosphere  duct. 
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2.   Integral  Representation  of  the  Fundamental  Solution 

The  Wiener-Hopf  Problem 

The  fundamental  solution  is  obtained  by  taking  the 
Fourier  transform  in  x  of  Eq.(l.7)  and  solving  the  resulting 
ordinary  differential  equation  in  y  which  formally  satisfies 
the  transformed  equivalents  of  Eqs.  (1.8  -  1.9)-    The 
transform  is  defined  for  convenience  to  be  decomposed  by  sum 
as 


G(y,a)  =  G^(y,a)  +  G_(y,a) 


(2.1) 


where 


+   00 


G^(y,a)  = 


+1 


(2-ir)' 


;e   dx 


(2.2) 


o 


Upon  taking  the  inverse  transform,  one  obtains  the 
representation 


&{^Q>yJx,y) 


{x^,y\x,y) 


'1^  o'^o 


+ 


(2ir)^ 


dae 


-lax 


-00 


cos 


hY(b-y)+ikZj  "^"'^T(b-y) 


coshyb  +  ikZ^  ^il^^^ 


(2.3) 


G(0,a) 


*  The  equivalent  condition  of  Eq.(l.6)  then  becomes 


lim   II  (y,a) 


y^+e 


t  (-C.) 


(2Tr)^ 
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where 


00 
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g^(x^,y^|x,y) 


2Tr 


\J 


da  e 


•la(x-Xo) 


coshyb+ikZ-p  ^— 


sinh  yy 
Y 


slnhyy^ 


cosh'y(b-y  )  +  IkZ. 


slnhY(b-yo) 


for  0  <  y  <  y^ 


coshY(b-y )+ikZ- 


sinhY(b-y ) 


for  y^  <  y  <  b, 


and  where  y  has  been  defined  by 


y 


(2.4) 


(2.5) 


The  Integrands  of  Eqs.  (2. 3-. 4)  are  even  functions  of  y 
and  hence  possess  no  branch  points. 

The  fundamental  solution  formally   satisfies  the 
boundary  conditions  given  by  Eqs.  (l.lO-.ll)  provided  the 
following  relation 


C(a)G^(0,a)  +  D(a)G_(0,a)  +  E(a)  «  0 


(2.6) 


holds  for  all  real  a,  where 


♦Once  the  properties  of  G(0,a)  are  determined  formally. 
It  will  be  Indicated  how  Eqs. ( 1. 7-. 11 )  are  satisfied 
rigorously. 
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C(a)  =  lk(Z^+Z_)cosh  yb  -i-  y  slnh  yb 


(2.7) 


D(a)  =  lk(Z(^+Z-p  )cosh  yb  +  y  slnh  yb 


^S^  I 


(2.8) 


and 


E(a)  = 


exp(lax^) 


(27r)2 


coshy(b-y^)+ikZj 


sinhy(b-y^) 


Y 


(2.9) 


and  where  quadratic  terms  In  the  Impedances  have  been 
neglected.  I.e.,  Z-Z  ^  Z-j-Z^  ^  0.   Provided  g  has  a  Fourier 
transform  and  gePl^  and  ge~P2^  are  absolute  Integrable  in 
the  respective  ranges  [0,+oo)  and  (-oo,0]  for  P-,,  p   positive 
and  real,  Eq.  (2.6)  is  analytic  in  the  strip 
-  min(P-,  ,  Im  k  ,  Im  h  )  <  Ima  <  mln(p  ,  Im  k  ,  Im  h  ), 
where  k  and  h  are  complex  constants  corresponding  to  the 
respective  zeros  of  C(a)  and  D(a)  which  lie  closest  to  the 
real  axis.   Further,  the  functions  G  (y,a)  and  G  (y,a)  are 
then  analytic  in  a  for  a  in  the  respective  domain 
Ima  >  -P-i   and  Ima  <  Pp. 

The  essence  of  the  Wiener-Hopf  method  is  to  decompose 
into  a  product  the  coefficients  of  G  (0,a)  and  G  (0,a)  in 
Eq. (2.6)  so  that 


C(a)  =  C_^(a)C_(a) 


(2.10) 


and 
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D(a)  =  D_^(a)D_(a)  (2.11) 


The  regions  of  analytlcity  of  C  (a)  and  C  (a)    functions 
are  Ima  >  -Imk  and  Ima  <  Imk  respectively,  and  those  of 
D  (a)  and  D  (a)  are  Ima  >  -Imh  and  Ima  <  Imh  respectively. 
A  second  decomposition  into  a  sum  is  performed  so  that 


\ 


mHjT  -  V°)  +  "-("'  (2-12' 


where  H  (a)  and  H  (a)  are  analytic  in  the  respective  domains 

Ima  >  -min(h  ,p  )  and  Ima  <  minfk  ,p  ),  where  p   is  that 
^  o   o  o   o  o 

singularity  of  E(a)  which  lies  closest  to  the  real  a-axis. 
Equation  (2.6)  may  then  be  written  as 

Cja)  D  (a) 

gi^  G^(0,a)  +  H^(a)  =  -  ^^   G_(0,a)  -  H_(a).     (2.13) 

The  functions  H  (a)  are  bounded  for  real  a  and  are 
0(a   )  as  Ia|  — ^  oo  and  the  ratios  C  (a)/D  (a)  and 
D  (a)/C  (a)  are  also  bounded  for  real  a  and  are  0(l)  as 
|a|  — ^00.   (Both  of  these  assertions  will  be  verified  in 
the  next  section.)   Further,  provided  the  above  assumptions 
on  g  hold,  the  Riemann-Lebesgue  lemma  implies  that  G  (0,a) 
vanish  in  a  strip  about  Ima  =  0  as  |a|  — ^  oo.   This  strip 
is  common  to  both  sides  of  Eq.  (2.13),  which  are  equal  to 
some  function  analytic  in  the  strip.   This  function  may  be 
continued  analytically  throughout  the  entire  a-plane,  but 
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because  of  the  above  described  asymptotic  behavior,  the 
resulting  entire  function  must  be  zero.   Therefore  one 
obtains  uniquely 

a(o,a)  =  -  If}  -  i.(Z,-Z3)  ^f^  'J^  (..») 

The  first  term  in  Eq.  (2.14)  is  0(a~  )  and  the  second  term 
Is  0(a~  )  as  |a|  — ^  oo  and  G(0,a)  is  bounded  and  continuous 
for  all  real  a,  provided  the  above  assertions  about 
C_(a)/D_(a)  and  H  (a)  can  be  demonstrated.   It  remains  to 
decompose  C(a),  D(a)  and  E(a)/C  (a)D  (a)  and  show  that  the 
assertions  may  be  verified. 


3.    Decomposition 

The  decomposition  into  a  product  given  formally  by 
Eqs.(2.10)  and  (2.II)  and  then  by  sum  gi\'en  by  Eq.(2.12) 
is  the  crucial  part  of  any  problem  solved  by  the  Wiener- 
Hopf  technique.   The  procedure  is  so  well-known  that  these 
factors  can  frequently  be  written  down  by  inspection.   In 
essence,  Eqs.  (2.7)  and  (2.8)  are  expressed  as  infinite 
products  and  Eq.  (2.12)  decomposed  by  an  application  of 
the  Cauchy  integral  theorem.   By  standard  methods  [see 
for  example  Knopp  (1957)]  one  obtains 
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C^(a) 


^o^ 


+ 


Yb 


1  -  i^f 


_  l_ab 


^*  iab/nir 


(3.1) 


n=o 


where  A  satisfies  the  equation 


A  IT   tan  A^TT  =  ikb(Z-r+Z-r  ) , 
n       n         I   L 


(3.2) 


and  by  successive  approximations  using  the  fact  that 
I Z^+Z, I  «  1,  one  obtains 


^o  =  i  [ikb(Zj+zJ]^ 


and 


A  =  n  +  — 
n        n 


n  =  1,2,3,. . 


(3.3) 


(3.4) 


The  other  coefficient  is  decomposed  as 


O  IT 

D  (a)  =  -^ 


o 


1  -  (^)'^ 


iab 


+  a^ir 


+  iab/mr 


(3.5) 


where   o   satisfies  Eq.  (3.2)  with  Z   replaced  by  Z  . 

Similar  expressions  to  Eq .  (3.3)  and  (3.'4)  may  all  be 

written  for  a  . 
n 

It  follows  from  Eqs.  (3-1)  and  (5.5)  that 

C  (a)        C_(a)   -fp 
a-^oo  +^^   a-^oo   -^^   n=l   n 


(3.6) 


Hence  C  (a)/D  (a)  and  D_(a)/C_(a)  are  0(l)  as  was  asserted. 
The  zeros  of  Eqs.  (3.1)  and  (3-5)  occur  at  ik^  and  +h^ 
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respectively,  where 

p      p      A  TT  2  Q      o      '^n'"'  2 

^n  =  ^   -  (-^)    ^"d    h2  =  k^  -  (^)^  (3.7) 


These  roots  lie  In  the  second  and  fourth  quadrants  of  the 

a-plane  and  off  the  real  axis  as  shown  In  Figure  2.   The 

condition  that  h   lies  closer  to  the  real  axis  than  k  Is 

o  o 


^l'  "  'S'  >  Im(fi^-fig) 


where 


fi^  =  Zj  +  Z^  +  ikb  (3.8) 


and 

fig  =  Zj  +  Zg  +  Ikb  (3.9) 


A  sufficient  condition  that  h  be  closer  to  the  real  a-axis 

o 

than  k  may  be  found  from  the  above  Inequality  to  be 


kb  +  Re(Z-j.+Z^)  >  I  fig  I 

which  Is  generally  satisfied  when  the  land  Impedance  Is 
much  greater  than  the  sea  and  Ionosphere  Impedance. 

It  has  already  been  mentioned  that  E(a)/D  (a)C  (a)  Is 
analytic  In  the  strip  -mln(lm  h  ,  Im  p  )  <  Im  a 
<  mln(lm  k  ,  Im  p  ) .   It  Is  therefore  possible  to  express 
this  function  as  the  Cauchy  Integral 
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E(a)        1   i,E(z)/D^(z)  C_(z) 

D^(a)C_(a)  "  2^  J  ^^^^ "^^  ^^'^^^ 

L 

where  L  Is  a  rectangular  path  shown  in  Pig.  5-   Provided 

Y     ^   0,    the  contributions  to  the  Integral  from  the  ends  of 

the  rectangle  vanish  exponentially  as  these  ends  are  taken 

to  Infinity.   Prom  the  remaining  portions  of  the  path,  the 

H  (a)  functions  are  obtained 

ic+oo 

E(z)/D  (z)  C  (z) 

^ dz  (3.11) 


«+^«)  =  2?r 


»-/ 


z-a 

IC-oo 


and 


H  (a)  =  -   ^ 


2Trl 


Id+oo 

^   E(z)/D,(z)  C  (z) 

+   dz  (3.12) 


z-a 
Id-oo 


where  -mln(lm  h.  Imp  )<c<Ima<d<  mln(lm  k  ,  Im  p  ) . 
To  evaluate  Eqs.  (3. 11-  .12)  the  line  Integrals  are 
displaced  Into  the  upper-half  z-plane  for  x  >  0  and  into 
the  lower-half  plane  for  x  <  0.   As  the  paths  are  displaced 
to  infinity  in  either  half-plane,  contributions  from  poles 
are  accrued  and  the  remaining  contribution  to  the  displaced 
path  becomes  exponentially  small.   Note  from  the  above 
that  the  H  (a)  functions  are  0(a   )  as  |a|— ^00,  as  was 
asserted  earlier. 
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In  what  follows,  x   is  taken  to  be  positive,  i.e.,  the 
transmitter  is  over  the  land  portion  of  the  duct.   Then 
upon  displacing  the  integration  path  into  the  upper-half 
z-plane  as  shown  in  Pig.  4,  one  obtains  for  H_(a) 


H_(a)  =  - -J  2_  Ir..}.    W   h   (i 


TTTTf  ^ —  (a-k  )k  b   ^k  +h  M   '^A/  k  +h, 
(2Tr)^  n=o  ^    n'^n     no. II-,    i        n      i 


+k  nn^  o.  k  - 

+h  M   ''A/  k  - 

n   o  JI-,    £  n 
^=1 


At      ikbZ      At 

^°s  ^-  ^o  +  a;;^—  ^^"  -^  ^o 


(3.15) 


where  quadratic  terms  in  the  impedances  have  been  neglected, 
When  a  4=  k  and  y     ^   0,    the  above  infinite  series  is 
uniformly  convergent.   The  H  (a)  function  is  then  found 
directly  by  means  of  Eq.  (2.12). 

Therefore  with  Eqs.  {^.l),    (3-5)  and  (3.13)  it  is 
possible  to  obtain  an  explicit  series  representation  for 
the  fundamental  solution  of  all  x  in  (-00,00)  and  all  y  in 
[0,b]  by  residue  calculus.   The  residues  are  due  to  poles 
at  ik  ,  ±h  and  ip  >n   =  0,1,2,3---,  where  P„  is  given  in 
terms  of  the  relation 

.M-^T"  .2 


P 


,2 

n    '^     "^  b 


k^  -  (-^)  (3.14) 


where,  in  turn,  the  [x     are  the  zeros  of  E(a)[    =  0  found 

"  ^o 

from  the  equation 
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fx^TT  cot  \i^ir   =  -IkbZ^  (3.15) 

which  has  approximate  solutions  for  |Z^1  «  1, 

^^n^i  +  —^.  (3.16) 

Using  arguments  similar  to  those  given  above,  it  can  be  shown 
that  generally  |  Im  h  |  <  |  Im  p  |  whenever  |Z<-,|  «  lZ-|. 

O  O  O  J. 

All  the  formal  results  are  now  available  to  allow  one 
to  show  that  Eq.  (2.3)  does  satisfy  the  conditions  given  by 
Eqs.  (I.7-.I3).   The  integrands  of  Eq.  (2.4)  and  that  given 
explicitly  in  Eq.  (2.3)  are  bounded  and  continuous  for  all 
real  a  and  vanish  exponentially  for  all  x  and  y  in  the 
respective  ranges  -00  <  x  <  00  and  0  <  y  <  b  except  for  a 
neighborhood  about  (x  ,y  )  in  the  former  case  and  a 
neighborhood  about  the  origin  (0,0)  in  the  latter  case. 
Then  exclusive  of  these  neighborhoods,  the  integrals  coverge 
uniformly,  limits  may  be  taken  under  the  integral  and  it  is 
routine  to  show  that  the  fundamental  solution  given  by  Eq. 
(2.3)  indeed  satisfies  Eqs.  (I.7-.II). 

In  the  neighborhood  of  the  shoreline  the  fundamental 
solution  is  given  by 

k(Z  -Z  )M 
g(xQ,y^|p  cos  e,p    sin  6)   =  i^^— p  cos  6    log  p 

+  a  bounded  function       (3-17) 
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where  M-,  Is  a  complex  constant  depending  upon  x  ,  y  and  the 
physical  parameters  of  the  problem.   Prom  Eq.  (3-17)  It  Is 
seen  that  the  condition  given  by  Eq.  (1.12)  is  satisfied. 

Equation  (3.17)  Is  obtained  by  separating  the 
Integration  Into  two  ranges:   a  range  bounded  away  from  a  =  0, 
from  which  the  logarithmic  term  Is  obtained,  and  the 
complement  of  this  range.   The  Integral  over  this  latter 
range  is  bounded.   This  coupling  between  the  behavior  at  the 
discontinuity  and  large  values  of  a  is  a  well-known 
phenomenon  (See  for  example  Noble  19?8).   Generally,  In 
diffraction  problems,  the  behavior  at  the  discontinuity  is 
specified  and  the  corresponding  condition  for  large  a  is  used 
to  separate  the  two  sides  of  Eq.  (2.13). 

For  Ixl  »  Ix  I  »  0,  the  modulus  of  the  dominant  term 
in  the  asymptotic  form  of  the  fundamental  solution  is  less 
than  Mp  exp[-| min(lmk  ,  Imp  )|x]  for  x  positive  and  less 
than  M^  exp[mln|(lmh  ,  Imp  )|x]  for  x  negative,  where  M^  and 
M^  are  bounded,  positive  and  real  constants.   Hence  it  can 
be  seen  that  Eq .  (l.lj)  is  satisfied.   In  a  similar  way, 
one  may  recover  the  fact  that  the  fundamental  solution  near 
(x  ,y  )  is  -( l/7r)log[  (x-x  )  +(y-y  )  ]^  plus  a  bounded  function 
which  is  the  standard  source  term  in  two-dimensional 
problems.   Derivatives  of  the  fundamental  solution  near 
(x  ,y  )  vary  as  the  inverse  distance  from  the  source  plus 
logarithmic  terms. 
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4.   Coastal  Refraction  of  VLF  Wave 

4. 1   Explicit  Series  Representation  of  Fundamental  Solution 
Let  J(y,a)  be  defined  as 


J(y,a)  =  cosh  r(b-y)  +  IkbZ^  3lnh^y(b-y) 


(4.1) 


Note  that  the  zeros  of  j(0,a)  are  those  of  E(a)j   ^  which 

'  -^o 
are  determined  from   Eqs.  (3.l4)  and  (3.I5).   Using  the 

above  equation,  Eq.  (2-3)  may  be  rewritten  as 


/     I    N    1   d  a  e  '"'"  '^o 
;(x^,yjx,y)  =  ^  jj^^-^ 


f     ^      -ia(x-xj  (  ^^^Y^  J^V«) 


O  "  o 


2Tr 


-00 


T 
sinh  yy 


for  0  <  y  <  y^ 


Y 


J(y.ct) 


v_ 


for  Yq  <  y  <  b 


00 


+ 


(2Tr)^ 


da  e J(y,a) 

J(0,a)c(a) 


laXr 


(27r): 


-  J(y^,a)  +  D  (a)C_(a)H_(a) 


+ 


00 

r 


(2t)2 


"(C)D(a^"^  [D^(a)C_(a)H_(a)] 


(^.2) 


The  above  integrals  may  be  evaluated  by  residue  calculus. 
All  integrands  in  Eq.  (4.2)  have  poles  at  ±p  .   The  Integrand 
of  the  second  integral  has  poles  at  JTk   and  that  of  the  third 
integral  has  poles  at  +h    (n  =  0,1,2...).   Thus  the  series 
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representation  for  Eq.  (4.3)  is 


g(xQ,yQ|x,y)  =  g2(xQ,yQ|x,y) 


+    (2t)2 


i  ^-  e      ^"  J(y,P^)H_(p^) 


n=o 


aj(0,a)/8a| 


n'    -^'-  n 


a=p 


n 


^(Pn) 


^-(Pn) 


+    (2Tr) 


00  -ihnX 

2 


^■^/^   J(y,hjD^(hJC_(h    )H_(h    ) 


n'    +'    n'    -'    n'    -'    n 


n=o 


J(0,  h    )D'(h^) 
\    '     n '       ^    n ' 


(4.4) 


for  X  <  0, 


g(xQ,yQ|x,y)   =   g2(x^,y^|x,y) 


+    (2t)2 


1  ^-  e        "  J(y,P^)H   (-p^) 


n  '    -  '    '^  n 


n=so 


Sj(0,a)/aa|      ^ 


n 


'^-(Pn) 


-     e-^^n(x-xo)   +  w/2   j(y,^^ ) j(y^,^J 

t 


n=o 


j(0,k    )C    (k    ) 


+    (2t) 


00 

1    -^: —   e 

2 


ik^x    +    l7r/2    J(y,k     )D    (kJC     (k     )H    (-kj 


n '    - ^    n '    +^    n ' 


n=o 


J(0,k    )c'(k    ) 
^    '    n '       ^    n ' 


('4. 5) 


for  0  <  X  <  x^  and 
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s(xQ,yQ|x,y)   -   g^{^^,Yj^,Y) 


+   (27r)^ 


^     .      ^lp,x   +   W2   j(y,pjH_(-p^) 


n=o 


SJ(0,a)/aa|^^ 


"n 


D(PJ 


n 


^ITpD 


n  ■ 


+ 


00        ikn(x-xo)    -r/      ,     vt-/        ,     \ 
e  J(y,k    )J(y   ,k    ) 


n=o 


J(0,k    )c'(k    ) 


+    (2t)^ 


i  ^  e'^"^  ^   '^/2   J(y,k    )D   (kJC    (kjH    (-kj 


n'    -^    n'   +^    n' 


n' 


n=o 


J(0,k    )c'(k    ) 


(4.6) 


for  X  >  X    ,    where 
o 


g2(x^,yjx,y)   = 


IPj^Ix-XqI    +   l7r/2 
dJ(0,a)/8a|       T" 


n 


J(y,Pn)J(yQ,Pn) 


sin   [i^ir  y/b 


M-n'"'/^ 


J^yo'Pn) 


+  < 


f or   0  <  y  <  y^ 


sm  H^TT  y/b 


M- 


n' 


1  for  y^  <  y  <  b 


(^.7) 


In  writing  Eqs.  (4. 4-. 7)  use  has  been  made  of  the  fact 
that  C(a),  D(a),  J(y,a)  and  j(y  ,a)  are  even  functions  of 
a  and  that  C  (a)  =  C_(-a)  and  D  (a)  =  D_(-a). 
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The  various  terms  In  the  series  representation  lend 
themselves  readily  to  physical  Interpretation.   The  first 
two  summations  In  Eqs.  (4. 4-. 6)  are  the  "sky-wave"  portion 
of  the  solution,  I.e.  that  part  which  Is  unobservable  If 
the  receiver  Is  on  the  ground  (y  =  O).   The  remaining 
summation  In  Eq.  (4.4)  are  "ground  vijaves"  transmitted 
beyond  the  shoreline.   The  last  two  summations  In  Eq .  (4.5) 
are,  respectively,  ground  waves  from  the  source  which  are 
Incident  on  the  shoreline  and  ground  waves  reflected  back 
from  the  coast. 

Waves  from  the  source,  but  receding  away  from  the 
shore  In  the  region  x  >  x  ,  and  the  ground  waves  reflected 
from  the  coast  are  given  respectively  as  the  last  two 
terms  of  Eq.  (4,6). 


4.2.    Reflection  and  Transmission  Coefficients 

When  the  receiver  and  the  transmitter  are  remote  from 
the  shoreline  and  each  other,  only  the  n  =  0  modes  of  Eqs. 
(4. 4-. 5)  are  observable.   To  this  remote  observer,  there 
exists  a  simple  wave  Incident  on,  reflected  from  and 
transmitted  beyond  the  discontinuity  at  the  shoreline. 
Weltz  and  Keller  (1955),  In  a  different  context,  have 
expressed  the  situation  where  a  simple  wave  propagates  past 
a  discontinuity  In  terms  of  reflection  and  transmission 
coefficients. 
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Consider  the  source  to  be  a  short,  two-dimensional 
dipole  at  (x-,  ,y-,  )  with  current  moment  IL  so  that  the 
current  density  is  given  as 


JyK'^o)  =  IL5(x^-x^)6(y^-y^) 


(4.8) 


A  vertically  oriented  antenna  at  the  receiver  observes  the 
E  component  of  the  electric  field.  For  the  above  source, 
this  field  component  is  given  as 


"^  o  1 


(^.9) 


For  I x-x  1  large  and  x  large  and  positive 


T-  /    i.\  -i£x     l|jc 

E  (x,y,t)  a^  e   ^  +  pe  ^ 


(4.10) 


and  for  x  large  and  negative 


E  •v  xe   ° 

y 


(4.11) 


where  f  is  either  p   or  k  ,  which  ever  has  the  smallest 
^  -^o     o 

imaginary  part.   It  is  assumed  that  min(lm  p  ,  Im  k  )  >  Im  h 
in  any  case.   The  reflection  coefficient  p  is  defined  by 


pTef lected 
P  ~  ^incident 


(4.12) 


x=»o 


and  the  transmission  coefficient  by 
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-.transmitted 


T  =-^ 


incident 


(^.13) 


x=o 


If  either  the  transmitter  or  the  receiver  are  near 
the  ground  then  only  the  mode  corresponding  to  ^  =  k 
Is  observable.   If  |  Im  pi  >  |  Im  k  |,  ,  then  ^  =  k  In 
any  case.   However,  when  | Im  k  |  >  | Im  p  |  and  y  and  y 
are  not  small,  |  =  p  .   The  condition  that  this  latter 
statement  be  true  Is  given  by 


Q. 


L' 


fi  I  >  Im(n  -fij) 


where  0,^    and  $7o  are  given  by  Eqs.  (3- 8)  and  (3-9) 
respectively  and 


Q,^   =  2Z^  +  Ikb  - 


ITT 


(h.lh) 


A  sufficient  condition  that  Im  k  >  Im  p   Is 

o  ^     -^o 


kb  +  Re(Z  +Zj)  >  IfijI 


which  Is  usually  satisfied  when  | Z, I  »  Iz,|  or  when 
I Z^  I  >  |Z.^|  and  the  earth-ionosphere  duct  is  a  quarter 
wavelength  in  height. 

Because  there  is  an  ambiguity  as  to  which  mode  Is 
propagating,  the  two  cases  are  discussed  separately.   For 
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case  I,  where  ^  =  k  characterizes  the  dominant  mode  for 

o 

X  >  0,  the  reflection  coefficient  p^  calculated  from 
Eqs.  (4.5)  and  (4.9)  Is  given  as 


D(k  )C  (k  ) 

where  H_(a)  =  -exp(lk  x  )H_(a)/(2Tr)^ . 

When  quadratic  terms  In  the  Impedances  are  neglected. 


p^  is  given  approximately  by 


h  -k 

Pi   h  +k  ^  -^o; 

o   o 


or  since  kb  »  |Z--+Z^|, 


Pl='^  •  ("-IT) 


The  transmission  coefficient  for  case  I,  calculated  from 
Eq.  (4. 4-. 5)  and  Eq.  (4.9)  is 

t         ^ 

h  D  (h  )C  (k  )c(h  )H  (h  )j(0,k  ) J(y, h  ) 

T-.  = 1 (4.10) 

k  C,(h  )D  (h  )j(y  ,k  )j(0,h  )j(y,k  ) 

0+^0'    ^  O   ^  "^  O   0^0   ^     o 

which  to  quadratic  terms  in  the  impedances  is  given  as 

2h. 


o 


I   h^+k 
o  o 


(4.19) 
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Equations  (4.16)  and  (4.19)  are  the  same  form  as  one 
classically  obtains  for  a  plane  wave  passing  between  two 
dissimilar  dielectric  media. 

For  case  II,  where  ^  =  p   characterizes  the  mode  for 
X  >  0,  the  situation  is  not  so  simple.   Using  Eqs.  (^-5) 
and  (4.9),  one  obtains  for  the  reflection  coefficient  Pjj, 


PlI  =  Pl^K)^l(Po'^o) 


(4.20) 


where  R(x  )  is  the  "scale  range"  defined  by 


R(x  )  =  exp[i(k  -p  )x  ] 


(4.21) 


and  P-,(p  ,k  )   is  of  polynomial  order  defined  by 


O  '  O' 


O   O '  +'  O 


sm^L^ry/b 


1  + 


-^  C(p  ) 


■^ r-  for  0  <  y  <  y 


Sinn  Ty/b       /   /  ^ 

j(y  ,p  )    f^°^  yo  <  y  <  ^ 

^  o   o ' 


-1 


(4.22) 

The  transmission  coefficient  for  case  II  is  found  from 
Eqs.  (4. 4-. 5)  and  Eq.  (4.9)  to  be 


T^^  =  1  +  p^R(xJpJp  .k  )F„(p  ,hj 


II 


o^  l^^o'  o^  2^^o'  O' 


(4.23) 
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where  P^fp  ,h  )  is  defined  by 


P^(p  ,h  )  = 


h^C(h^)D^(h^)H_(h^)j(y,h^)aj(0,a)/Sa 


p  J(0,h  )C  (h  )H  (-P  )j(y,p  )D  (h 


a  =  P^ 


D  (p  )    C  (p  ) 


C  (p  )  "  D    (p  ) 
-   O  '      -f^  O  ■ 

(^.2^) 


However,  from  Eq.  (4.21)  it  is  seen  that  when  Im  k  >  Im  p 
the  scale  range  vanishes  when  |x  |  — ^00.   Therefore  it  is 
possible  to  conclude  that  the  dominant  mode  is  totally 
transmitted  and  no  reflected  component  is  observable. 
Even  though  p   characterizes  the  dominant  mode  in  this  case, 
any  reflection  must  come  from  the  k  mode.   Provided  the 
impedance  of  the  ionosphere  is  sufficiently  small  and  the 
source  and  the  receiver  are  elevated  and  far  from  the 
shoreline  and  each  other,  the  discontinuity  is  not  . 
observable  to  lowest  order  in  the  impedances.   On  the 
other  hand,  if  the  ionosphere  impedance  is  not  small,  or 
if  either  the  source  or  receiver  is  near  the  earth's 
surface,  the  portion  of  the  incident  wave  reflected  is 
proportional  to  the  difference  of  the  land  and  sea 
Impedances. 
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4.5    The  Recovery  Effect 

Strictly  speaking,  because  of  the  complex  Impedances, 
every  mode  described  by  Eqs.  (4. 4-. 7)  Is  propagating, 
though  those  of  higher  order  (n  >  O)  are  strongly  damped 
at  low  frequencies.   At  slightly  higher  frequencies, 
however,  more  than  one  mode  may  be  sustained,  particularly 
over  the  sea.   At  these  higher  frequencies,  an  electric  field 
consisting  of  one  dominant  mode  and  an  Infinite  number  of 
heavily  damped  modes  is  incident  upon  the  shoreline;  the 
field  transmitted  beyond  may  consist  of  more  than  one 
observable  mode.   The  energy  from  the  supressed  modes 
over  land  would  have  been  transferred  to  the  higher  order 
modes  over  the  sea.   This,  in  essence,  is  the  "recovery 
effect"  for  VLP  transmission. 

Mllllngton  (1949)  noted  that  the  recovery  effect  was 
strongest  when  the  observer  was  near  the  earth's  surface. 
The  same  conclusion  may  be  Inferred  from  Section  4.2 
where  It  was  noted  that  reflection  from  the  shoreline 
discontinuity  was  observed  when  the  observer  or  the 
transmitter  was  near  the  earth.   Therefore  in  order  to 
Illustrate  the  recovery  effect,  the  observer  is  assumed 
to  lie  in  a  region  near  the  earth  and  in  a  range 
(im  k-,  )~  <  |x|  <  (Im  h-,  )   .   Such  a  region  is  shown  as 
shaded  in  Pig.  5.   In  this  case,  only  one  mode  Is 
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observable  for  x  >  0  and  two  are  observable  for  x  <  0.   Over 
the  land  (x  >  O),  the  y-component  of  the  electric  field 
behaves  like  Eq.  (4.10)  while  over  the  sea  (x  <  O)  one  has 

r    -Ih^x      ,         V  -ih-ix,        ,,    . 
Eya^Tj[e   °  +  P^(h^,h^,k^)e    ^  ]       (^.25) 

for  |x|  in  the  prescribed  range.   The  function  P^(h  ,h,,k  ) 
is  given  as 

h  J(y,h  )c(h  )D  (h  )H  (h  ) 

P^  h^,h,  ,kj  =  -i i ^   +   ^  .V    ^  C^-26 

^        ^     h  J(y,,h^1c(h  )D  (h  )H  (h  ) 

O  ^^1'  O'  ^  o'  +^  o'  -^  o 

j(0,h^)D'(h^)C_^(h^) 
j(0,h^)D'(h^)C^(h^) 

Equation  (4.26)  may  be  written  approximately  as 


P  (h^,h^,k^)  ^  -ikb(Z^+Z^)  ^  ^  (4.27) 

oil 


where  quadratic  terms  in  the  impedances  have  been  suppressed. 

In  order  to  compare  the  field  on  either  side  of  the 
shoreline  for  the  range  |x|  >  (im  k-,)"   the  squares  of  the 
amplitudes  of  the  two  cases  are  calculated  and  found  to  be 

|Ey(x>0)l2  =  e''"^°\  |pj|2e-2^"^°^ 

(4.28) 

+  2|  pjl  cos(2Rek  x+argp-j.) 
for  an  over-land  observer  and 
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/     ^m2    I     i2)  2Imhox   ,   ,         m2  2Iinh-,x 
E  (x  <  0)r  =  |1+PtI   1  e     °  +  IP  (h  ,h  ,ki)|  e     1 


yK^^   ^    -/,      i-.t^j,   .  -        -  ,.^,..^,..^ 


+  2|F.(h  ,h,,kJ!e^'^^'^°+^l^''  cos[Re(h  -hjx 


+  arg  P^(h^,h^,k^)] | 


('!.29) 

for  an  over-sea  observer.   By  comparing  Eqs.  ('4.28)  and 

(4.29)  it  can  be  seen  from  Eq.  (4.2?)  that  the  effective 

Increase  in  amplitude  as  the  observer  moves  across  the 

shore  line  Is  proportional  to  |Z-p+Z^|  .   The  recovery  of 

1   J-j 

2 
E    Is  Illustrated  In  Pig.  5.   It  should  be  noted  that 

I    I  2 
while   E    Is  increasing  as  the  observer  moves  from  land 

I     yl  6 

I   I  2 
to  sea,  |E  I   Is  decreasing  by  an  amount  proport.i  ona  I  to 

I  Z-r-+Z„  I /I  Z-|-+Z-|.  I  .   As  a  result,  the  wave  vector  may  be 

thought  of  as  "straightening  up"  as  the  wave  passes  from 

the  relatively  poor  conducting  land  to  the  more  highly 

conducting  sea. 
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Figure  Captions 

1.  The  earth-ionosphere  duct  showing  schematic  placement 
of  transmitter  T  and  receiver  R. 

2.  The  alpha -plane  showing  the  occurrence  of  the  +k  terms. 

3-   Path  of  Integration  for  decomposition  by  sum  of  the 
E(a)/D  (a)C_(a)  function 

4.   Evolution  of  H  (a)  by  displacement  of  Integral  path  to 
Infinity  In  the  upper-half  z-plane. 

5-   Region  of  observation  of  the  recovery  effect  so  that 

one  mode  Is  observable  for  x  >  0  but  two  are  observable 

for  X  <  0. 

2 
6.   The  "recovery"  of  | E  |   as  the  observer  passes  from 

land  to  sea. 


-  32  - 


5.0    o 


IONOSPHERE 


■> 


R  • 


•T 


SEA 


LAND 


FIGURE   1 


55  - 


Im  a 


"o* 


K 


a- PLANE 


F(z.+z,) 


Re  a 


FIGURE    2 


-   54  - 


Imz 


STRIP  OF  ANALYTICITY 


PATH  OF  INTEGRATION 


r-PLANE 


•a 


min(Imko,Impo) 


-d 


Rit 


-c 


—  mln(Imho,Impo) 


FIGURE  3 


-   55  - 


Imz 


ORIGINAL  PATH 


z- PLANE 


O'sPLAceo  p^rn  of 


"^^^GRAT^ 


OF  INTEGRATION 


Rez 


FIGURE  4 


-  36  - 


(Imkj)"*  (Imhif* 


FIGURE  5 


-   37   - 


U+F3(hpho)l 


FIGURE    6 


-    38    - 
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such  contractor. 
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